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t:j- . Abstract 

> 

ly-v ' The general relativistic gravitomagnetic clock effect consists in the fact that two massive test bodies 

^D ' orbiting a central spinning mass in its equatorial plane along two identical circular trajectories, but in 

opposite directions, take different times in describing a full revolution with respect to an asymptotically 



<^ ', inertial observer. In the field of the Earth such time shift amounts to 10 ^ s. Detecting it by means 

^—^ ' of a space based mission with artificial satellites is a very demanding task because there are severe 



constraints on the precision with which the radial and azimuthal positions of a satellite must be known: 
6r ~ 10~^ cm and 6ip ~ 10^^ milliarcseconds per revolution. In this paper we assess if the systematic 



bJO' errors induced by various non-gravitational perturbations allow to meet such stringent requirements. 

A couple of identical, passive laser-ranged satellites of LAGEOS type with their spins aligned with 
the Earth's one is considered. It turns out that all the non vanishing non-gravitational perturbations 



^H ' induce systematic errors in r and if within the required constraints for a reasonable assumption of the 

mismodeling in some satellite's and Earth's parameters and/or by using dense satellites with small 
area-to-mass ratio. However, the error in the Earth's GM is by far the largest source of uncertainty 
in the azimuthal location which is affected at a level of 1.2 milliarcseconds per revolution. 



1 Introduction 

General Relativity, in its weak-field and slow-motion approximation, predicts that two massive 
test particles following in opposite directions two geodesic circular orbits of identical radius 
in the equatorial plane of a central spinning body take different times in describing a full 
revolution with respect to an asymptotically inertial observer. If t"*" denotes the period of the 
particle moving in the same sense of the rotation of the central body and t~ denotes the period 
of the particle moving in the opposite sense, it turns out that 

At = t+-t-=47r^, (1) 

where J and M are the proper angular momentum and the mass, respectively, of the central 
source and c is the speed of light in vacuum. Note that the Newtonian gravitational constant 
G does not appear in eq.(|l|). The co-rotating particle turns out to be slower than the counter- 
rotating one. This is the so called gravitomagnetic clock effect [Cohen and Mashhoon, 1993; 
Mashhoon et al., 1999; 2000]. Eq.(|T]) is also equal to the difference in the test particles' proper 
periods r^ for test particles orbiting at distances larger than the gravitational radius r^ = 
1GMj(? of the central body and neglecting terms of order 0[c~'^\ 

For a couple of artificial Earth satellites eq.(|l|) yields a time shift of 10"^ s which is large 
enough to be detected according to the present-day level of accuracy in timing measurements 
with atomic clocks [Lichtenegger et al, 2000]. Note that eq.([l|) is independent of the satel- 
lite's orbit radius. This feature could be fruitfully exploited in choosing suitably the orbital 
parameters of the satellites to be employed. However, it has been shown in [Gronwald et al, 
1997; Lichtenegger et al, 2000] that, in order to make feasible such measurement, the radial 
and azimuthal positions of the satellites should be known at a level of 6r ~ 10~^ cm and 
6ip ~ 10^^ milliarcseconds (mas) per revolution. Such constraints are very stringent: suffice it 
to say that the position of LAGEOS laser-ranged satellite is presently known at a centimeter 
level. Recently, in [lorio, 2001; Lichtenegger et al, 2001] the impact of the long-periodic and 
short-periodic gravitational perturbations, with the systematic errors induced by them in r and 
if, have been worked out: the conclusion is that, at the present level of knowledge of the static 
and time- varying parts of the Earth's gravitational field [Lemoine et al, 1998], the gravitational 
errors are larger than the required 6r and 6ip. 



However, in the near future the new, more accurate data for the terrestrial gravitational 
field from the CHAMP and GRACE missions will be available and the situation could become 
more favorable. So, it appears important to assess the error budget due to the non-gravitational 
forces [Milani et al., 1987] on the satellite's radial and azimuthal locations, which is the scope 
of this paper. 

We will consider throughout it a couple of identical passive, geodetic, laser-ranged spherical 
satellites of LAGEOS type both because Satellite Laser Ranging (SLR) has reached in the 
last decade an astonishingly level of accuracy and because with this kind of satellites it is 
by far simpler to model accurately enough the non-gravitational perturbations. Indeed, they 
depend on the physical and geometrical characteristics of the satellites and on the geometry 
of their orbits as well. For example, twice a year, almost six months apart, the Sun, moving 
along the ecliptic, intersects the Earth's equatorial plane so that the satellite's revolutions 
occurring these times are affected by the phenomenon of eclipses. The radiative perturbations 
which are generated, directly or indirectly, by the solar electromagnetic radiation like the direct 
solar radiation pressure, the albedo and the Yarkovsky-Schach effect act differently on the 
satellite's orbit with respect to the orbital arcs described in full sunlight and in many cases 
such discrepancies are relevant and difficult to calculate. Moreover, the thermal forces which 
cause the Yarkovsky-Schach and the Rubincam effects depend on both the physical properties 
of the satellite and on the orientation of its spin axis with respect to an inertial frame. The 
terrestrial environment may cause the spin direction to change in a more or less predictable 
way over time spans of some years, as it is occurring to LAGEOS [Metris et al., 1999]. 

In order to detect the gravitomagnetic time shift, which is cumulative, there is no need of 
very long orbital arcs: indeed, for example, the orbital period of LAGEOS amounts to 1 x 10~^ 
days only, so that 10^ — 10'^ revolutions [Lichtenegger et al., 2000] would correspond to arcs 
10^ — 10^ days longQ. We could take advantage of this fact by choosing the arcs so to avoid 
the eclipses effects; moreover, over such short time spans it would be possible to consider the 
satellite's spin axis direction as fixed in the inertial space. Indeed, many thermal effects vanish 
for a suitable fixed spin axis direction. 



^However, as pointed out in [lorio, 2001], short arcs would not allow to average out many gravitational 
perturbations 



The basic assumptions of our study are the following: 

• In order to simplify the calculations we will consider a couple of LAGEOS type satellites with 
semimajor axis a = 12270 km, as for LAGEOS, zero eccentricity e and inclination i and orbital 
period P = 1.35 x lO'' s 

• The physical parameters of the satellites like area-to-mass ratio S/m, reflectivity and drag 
coefficients Cr and Co, etc. are assumed to be equal to those of LAGEOS: S/m = 6.87 x 10"^ 
cm2g-i, Cfi= 1.13, Czp = 4.9 

• We will assume that their spins are aligned with the Earth's spin axis assumed as z axis of a 
terrestrial equatorial inertial frame 

• Only orbital arcs in full sunlight will be considered: the eclipses effects will be neglected 

• The perturbations on r and ip will be averaged over an orbital revolution 

• The perturbations which will be considered are the direct solar radiation pressure, the Earth's 
albedo, the Poynting-Robertson effect, the direct Earth IR radiation pressure, the anisotropic 
thermal radiation or Yarkovsky-Schach effect, the thermal trust or Rubincam effect, the atmo- 
spheric drag and the effect of the Earth's magnetic field. 

The paper is organized as follows: in Section 2 and Section 3 the perturbations on r and ip, 
respectively, are worked out, while Section 4 is devoted to the conclusions. 

2 The radial position 

The radial position of a satellite in a perturbed circular orbit may change due to variations 
both in its semimajor axis a and its eccentricity e according to [Christodoulidis et ai, 1988]: 



Ar = W(Aa)2 + -{aAey. (2) 

The perturbations on a are calculated with [Milani et ai, 1987]: 

a = -T, (3) 

n 

where n = ^ is the mean motion, P is the Keplerian satellite's orbital period and T is the 
along-track disturbing acceleration. Concerning the eccentricity, since we are dealing with 
circular orbits we will use the components of the eccentricity vector h = ecoscu, k = esinu; {u 
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is the argument of perigee) and [Milani et al, 1987]: 

h = ecosto + 0(e) = —[-Rcos{u + f) + 2T sm{uj + f)] + 0{e), (4) 

na 

k = ecosuj + 0{e) = —[Rsm{uj + f) + 2T cos{uj + /)] + C»(e), (5) 

where / is the true anomaly and R is the radial component of the perturbing acceleration. For 
a circular orbit, from eq.(|D and eq.@ it follows: 

e = sJ{hY + {kY. (6) 

2.1 The radiative perturbations 

The perturbing acceleration due to direct solar radiation pressure can be approximately written 
as [Lucchesi, 1998]: 

Wq = -WqS, (7) 

where s is the unit vector from the Earth to the Sun and: 

We = --Cn. (8) 

m c 

In eq.(^ the term ((^0/^0)^, in which d^, is the semimajor axis of the Earth orbit around the 

Sun and Tq is its instantaneous distance, has been set equal to one. Iq is the solar constant. 

For LAGEOS type satellites eq.(H) amounts to almost 3.6 x 10^'^ cm s^^. 

Concerning the semimajor axis, it turns out that if the total solar radiation force acting on 
the spacecraft can be expressed in the general form F(s), no long-term effect in a will appear 
to any order in e [Milani et al, 1987]. So, < a >2-k= 0. 

Regarding the eccentricity, it can be proved that it is affected, at zero order in e, by long- 
term periodic perturbations with an almost yearly period [Lucchesi, 1998]. The amplitudes for 
the components of the eccentricity vector are: 

3 W0 , , 

2 na 

< k >2.oc 1^, (10) 

/ na 

so that, over an orbital revolution, the radial position of the satellite changes of almost 11 cm 
[Milani et al, 1987]. 



Assuming that the solar constant is known at 0.3% [Ciufolini et ai, 1997], the major source 
of uncertainty in the satelhte's radial position due to the direct solar radiation pressure would 
reside in the Cr satellite's coefficient. Assuming a global 0.5% mismodeling, the systematic 
error would amount to Stsrp = 5 x 10~^ cm. This result suggests that, in order to meet the 
requirement of 6r ~ 10^^ cm, a careful analysis of the optical properties of the satellites should 
be conducted in accurate pre-launch controls. 

For a circular and equatorial orbit the perturbations on the satellite's radial position induced 
by the Earth's albedo are entirely due to the changes in the eccentricity vector as well. Indeed, 
the along-track acceleration, which is particularly sensitive to the specular part of Earth's 
albedo and to its spatial and temporal variations, can be expanded in terms of long-periodical 
harmonics whose coefficients are proportional to sini [Anselmo et ai, 1983]. Concerning the 
eccentricity vector, the simple analytic model by Metris et ai, [1997], based on an uniform 
mean albedo of ^ = 0.3, turns out to be adequate. The related perturbation amounts to 
8 X 10^^ cm per orbit. Then, the systematic error due to the mismodeling in A, assumed to be 
10% [Lucchesi, 1998], amounts to 8 x 10^'^ cm per orbit. 

It may be important to note that a way to reduce the impact of the direct solar radiation 
pressure and the albedo could be the use of a couple of dense satellites with small — . 

The Poynting- Robertson acceleration [Burns et ai, 1979] leaves unaffected the eccentricity 
while changes the semimajor axis with long-term perturbations whose nominal amplitudes are 
of the order of 10^^ — 10^^ cm per orbit. Consequently, the related systematical errors are 
negligible. 

2.2 The thermal perturbations 

Concerning the direct Earth IR radiation pressure, the Earth thermal emissivity S can be 
modelled in a form of a latitude-dependent spherical harmonic expansion on a spherical Earth 
surface [Sehnal, 1981]: 



S = So + S,{t) Pi(sin0)+^2i'2(sin0) + ..., (11) 



where is the terrestrial latitude. Regarding the effects of the ffist two zonal constant terms So 
and £2, < a >2n= since it is proportional to (esini)^ and < e >2tt= since it is proportional 
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to e(sini)^ [Sehnal, 1981]. The perturbations due to Si{t), which shows an approximately yearly 
variation, vanish as well since R and T are proportional to sini [Metris et al, 1997]. 

The Yarkovsky-Schach effect [Afonso et al., 1989] induces perturbations on the semimajor 
axis which vanishes since the first non-zero term in a is of order 0{e) for orbital arcs in full 
sunlight. As far as the eccentricity vector is concerned, e would be perturbed at a level of 
9 X 10"^ cm per orbit. However, such effect depends on the orientation of the satellite's spin 
axis C,'- for a fixed direction such that ^x = ^y = 0, as it could be feasible over short arcs, 
< h >2tt=< k >2tt= [Lucchesi, 1998]. 

The thermal thrust or Rubincam effect [Rubincam, 1987] depends on ^ as well; for ^.j. = 
^y = it vanishes because it can be proved that, in this case, R and T are proportional to sini. 

2.3 Other perturbations 

The neutral and charged drag has negligible effects on r. Indeed, by neglecting the corotation 
of the exosphere, it turns out that [Milani et al, 1987]: 

S 
< Aa >2n= -Co—a^Q, (12) 

m 

in which q is the atmospheric density. For q = 8.4 x 10^^^ g cm^^ [Afonso et al, 1985] we 
would have a decay in a of 4.2 x 10"'^ cm per revolution. The same holds for the charged drag 
[Rubincam, 1990]. For circular orbits the eccentricity is not affected by the drag [Milani et al, 
1987]. 

If the satellites carry an electric charge, as it is the case for LAGEOS, the Earth's magnetic 
field acts upon them via the Lorentz force. Its effect on the semimajor axis is of course zero since 
the Lorentz force does not change the satellite's total mechanical energy W and, consequently, a: 
indeed, W = —^^- Regarding the eccentricity vector, it turns out that < h >2n=< k >2-k= 
because R is constant and T is of order Oie). 

3 The azimuthal position 

As in [lorio, 2001; Lichtenegger et al, 2001], for a satellite in an equatorial and circular orbit 
the rate of the azimuthal angle can be calculated by means of: 

(f = iu + tlcosi + M. (13) 



In eq.(|T3|) il is the longitude of the ascending node and 7V1 is the mean anomaly whose rate 



equation is given by: 

2 r 



M=n R--Vl-e^{dj + cosin). (14) 

na a 

For e = i = eqs.(|13|)-([l|) yield: 

^ = n-—R. (15) 

na 



In dealing with eq.([l5|) the indirect effects on n induced by the perturbations on a must be 
considered as well. At the first perturbative order they are: [Milani et al., 1987] 

An = -— Aa. (16) 

3.1 The radiative perturbations 

The direct solar radiation pressure does not affect if because, as seen in previous section, 
< CL >2tt= 0; moreover, < R >2-k= 0. 

Regarding the Earth's albedo, the indirect perturbations on n vanish because < a >2-k= 0. 
According to the model by Metris et al, [1997] < R >2-n= 0. 

The non- vanishing perturbations of order 0{e^) induced by the Poynting- Robertson effect, 
which amount to almost 10~^ mas per revolution or less are negligible. 

3.2 The thermal perturbations 

Concerning the direct Earth IR radiation pressure, < a >27r= so that the indirect perturba- 
tions on the mean motion vanish as well. The first two constant zonal terms of the Earth IR 
emissivity yield non- vanishing terms of zero order in e. Indeed, it turns out that, for £q: 

<--R>?}=--{^fSo-^. (IT) 

na na a m c 

which yields —1.6 x 10^^ mas per revolution. The systematic error induced by the mismodeling 
in £q is within the limit 5ip ~ 10^^ mas per revolution and could be reduced by using satellites 
with small area-to-mass ratio. The bias due to £2 is negligible since its nominal perturbation 
amounts to 1 x 10^"^ mas per revolution. As in for r, also in this case £\{t) does not contribute 
since R is proportional to sini. 
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The Yarkovsky-Schach effect does not affect the azimuthal position since < a >2-k= and 
it turns out that R averages out over an orbital revolution. The same holds for the Rubincam 
effect which is not present when ^^ = ±1 and i = 0. 

3.3 Other perturbations 

The indirect perturbation on n due to the drag shift experienced by the semimajor axis is 
negligible because it amounts to 6 x 10^^ mas per revolution, while: 

2 S 

< R>2t,= Co—nae KsixiE >2n=^-, (18) 

na m 

where E is the eccentric anomaly. 

The effect of the Earth's magnetic field is completely negligible since it is of the order of 
10^^ mas per revolution. 

The major source of error in the azimuthal position turns out to be the Earth's GM. 
Indeed, Acp^"^^^^ = ^^^3 ' x P = 1.2 mas per revolution by assuming 5{GM(f^) = 8 x lO"*^^ cm^ 
s-2 [McCarthy, 1996]. 

4 Conclusions 

In the context of the gravitomagnetic clock effect, by considering a couple of identical SLR 
satellites of LAGEOS type with e = i = and ^x = ^y = 0, ^z = =tl over short orbital arcs in 
full sunlight, it turns out that some non-gravitational perturbations on r and (/? vanish. 

Regarding the radial position, the largest perturbations are due to the direct solar radiation 
pressure which would change the satellite's distance of almost 10^ cm per revolution and the 
Earth's albedo which would induce a change of 8 x 10~^ cm per revolution. However, their 
systematic errors induced by the mismodeling in the optical properties of the satellites and 
the Earth's albedo should fall below the cutoff 6r ~ 10~^ cm. The influence of the other 
non vanishing perturbations, like Poynting-Robertson effect and neutral and charged drag, is 
negligible. The thermal perturbations vanish. 

The azimuthal angle is perturbed by the Earth's IR direct radiation pressure at a level of 
10^^ mas per revolution. The other non vanishing perturbations are negligible. 



It should be pointed out that all the non- vanishing perturbations are proportional to — . 
This means that the impact of their systematic errors could be reduced by using particularly 
dense satellites with small area-to-mass ratio. However, the largest source of error in ip is 
the uncertainty in the Earth's GM which induce a bias of 1.2 mas per revolution. This is 
a hard limitation to overcome because it is independent of the particular satellite employed 
and is related to our knowledge of the terrestrial gravitational field. Moreover, as pointed 
out in [lorio, 2001], at the present level of knowledge of the terrestrial gravitational field, the 
mismodeled gravitational time varying perturbations induce systematic errors which are larger 
than the required constraints; they could be overcome by averaging them over adequately long 
time spans. 

The results of the present work suggest that for a suitable choice of the satellites to be 
employed and of their orbital geometry it would be possible to keep the non-gravitational 
perturbations within the required constraints in order to make feasible the measurement of 
the gravitomagnetic clock effect. As far as the systematic errors induced by the forces acting 
upon the satellites are concerned, the major problems come from the Earth's gravitational 
environment. Improvements in satellite tracking accuracy of almost two orders of magnitude 
are needed as well. 
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